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Abstract 

Extending phase space to include time and it canonical conjugate 
energy as well as the usual momentum and position variables, and then 
r~| ! introducing the constraint which sets energy equal to the Hamiltonian, 

gives a symplectic action of the additive group R which corresponds to 
time translation along the solutions to the equations of motion. This 
allows the BRST operator of a system with secondary constraints to 
be constructed from first principles. 



1 Introduction 

This paper derives from first principles the extra ghost terms which must 
be added to the classical Hamiltonian in the BRST quantization of systems 
with secondary (and higher) constraints. 

BRST cohomology provides a systematic approach to the path integral quan- 
tization of systems with symmetry, effectively reducing the integrals to the 
true phase space of the system even when this is a finite-dimensional mod- 
uli space. In the Hamiltonian setting the symmetries of a system manifest 
themselves as constraints on the phase space, as was first shown by Dirac 
PP. For systems with first class constraints all of which are primary, the 
BRST method is well understood and can be justified in terms of Hamilto- 
nian group actions on symplectic manifolds [2J [3J HI [5] . However for many 
realistic systems the situation is more complicated, in particular when a sys- 
tem has secondary constraints not only must the correct BRST operator be 
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constructed, but also an extension of the Hamiltonian is required so that it 
commutes with the BRST operator. This paper extends the understanding 
of BRST quantization to systems with secondary (and higher) constraints, 
deriving the extra terms in the Hamiltonian from first principles (as opposed 
to simply adding terms so that the Hamiltonian commutes with the BRST 
operator) by extending phase space to include time and its canonical conju- 
gate energy. 

In the Lagrangian setting it was recognised by Noether that a conserved 
current could be associated with each symmetry generator of a system. A 
description of these ideas and their historical development may be found in 
Kosmann-Shwarzbach [6]. In his pioneering work on the canonical quantiza- 
tion of systems with symmetries [1] Dirac recognised that in the canonical 
setting symmetries manifest themselves as constraints on phase space. In 
the simplest cases a Lagrangian system with n degrees of freedom whose 
symmetries give rise to m Noether currents leads to a phase space with m 
primary constraints which are both first class and conserved in time; these 
primary constraints correspond directly to the Noether currents, and lead via 
the Marsden-Weinstein reduction process [7] to the true phase space of the 
system which has dimension In — 1m. However it was observed by Dirac that 
systems occur whose primary constraints are not constants of the motion and 
thus do not directly correspond to the Noether currents of the theory. The 
time variation of the primary constraints are normally introduced as sec- 
ondary constraints, but at the expense of the standard machinery which has 
been developed using constraints to understand the reduced or true phase 
space of the theory and hence proceed to quantization. 

This paper provides an understanding of such systems, by considering the 
space of phase space trajectories rather than simply phase space, building 
on the work of Castellani, Henneaux, Tetelboim, Zanelli, and Pons [H [9j HO] ■ 
A key point is that it is not appropriate to regard primary and secondary 
constraints as individual constraints on phase space, instead each primary 
constraint together with any corresponding secondary constraint must be 
built into a single object on the space of trajectories in phase space; these 
composite objects correspond to the momentum map of a Hamiltonian group 
action on the space of phase space trajectories. In more mathematical lan- 
guage, the theory is developed in terms of the space of paths in a symplectic 
manifold; this space itself has a symplectic structure, and it will be shown 
that more general gauge symmetries, of the kind which lead to secondary 
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(and possibly tertiary and so on) constraints can be formulated by extending 
standard constructions in Hamiltonian dynamics to the symplectic manifold 
of paths. 

A further important step is the extension of phase space to include time and 
its canonical conjugate energy, as described in Woodhouse [TT]. Imposing the 
constraint that the energy is equal to the Hamiltonian function then allows 
us to take a 'snapshot' in the space of classical trajectories in a manner which 
respects the dynamics of the system, so that a finite-dimensional phase space 
can be used rather than the infinite-dimensional space of trajectories, as is 
explained in section |5j 

In section [5] the standard BRST construction for systems with primary first 
class constraints is briefly reviewed. Important steps towards understanding 
the structures which arise when secondary constraints occur made by a num- 
ber of authors including Castellani, Henneaux, Tetelboim, Zanelli, and Pons 
[SJ EJ ED] are is described section [3] with particular emphasis on the work of 
Pons who recognised that Dirac had restricted his analysis to an infinitesi- 
mal neighbourhood of the initial conditions, while it was in fact necessary to 
consider the gauge parameters and phase space trajectories for all values of 
the time parameter t. 

In section HJ using rather more mathematical language and largely following 
[12], we review the symplectic geometry of the space of paths in phase space, 
and show how a Hamiltonian group action on this space leads to a momentum 
maps that encodes primary and secondary constraints in a single object and 
allows a finite-dimensional formalism to be used. In section the BRST 
operator, as well as the extended Hamiltonian, is constructed for systems 
with secondary constraints, and the full quantization procedure described. 
This is achieved by extending phase space to include time t and energy E, 
together with the additional constraint H — E = which corresponds to time 
translation along the solutions to the equations of motion. 

To set the scene we will describe a toy example introduced in [13] which gives 
a simple illustration of the issues involved and establishes notation. (Several 
more examples are analysed in the final section.) Consider a system with 
two degrees of freedom g l , i — 1, . . . , 2 and Lagrangian 




(1) 
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The action is invariant under transformations 



5{q 2 ) = e and = -e (2) 

for arbitrary e(t) and the equations of motion are 

q 2 + q 1 = and q 2 + q 1 = . (3) 

The Noether current corresponding to the invariance (T5]) is 

J(q i ,q i )=q 2 + q 1 , (4) 

which is readily seen to be conserved provided that the equations of motion 
hold. 

Passing, via the Legendre transformation, to the Hamiltonian formulation of 
this model we obtain 



9 £ n <9£ , 

Pl = W = P2 = W = q+q ' 



so that we have one primary constraint (f)^ = pi — 0. The Hamiltonian of 
the system is 

H(p^q l ) = +\p 2 2 -P2q l (6) 
from which we see that the equation of motion for the constraint (f)^ = p\ is 



Pi = V2 , (7) 

so that this constraint is not preserved in time unless the secondary constraint 
cj)^ = p\ = p 2 = also holds. This secondary constraint is preserved in time 
so that there is no tertiary or higher constraint. One purpose of this paper 
will be to reconcile the apparent mismatch between the number of Noether 
currents (in this case one) and constraints (in this case two) by building 
primary and secondary (and where they exist tertiary and so on) constraints 
into a single object. 

Throughout this paper we assume that there are no second class constraints 
or reducible symmetries. To include them would introduce further complica- 
tion without throwing more light on the developments in this paper. We also 
use the language of quantum mechanics rather than field theory, although 
field theoretic examples are considered in section |6j 
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2 Review of case where there are no secondary 
constraints 



For systems which simply have primary first class constraints there is a full 
understanding of the role of symmetry in the classical dynamics which allows 
path integral quantization of the system using BRST techniques. These ideas 
are reviewed in this section, using the formalism of quantum mechanics. 

Suppose that an n- dimensional system with position variables q l {t),i = 
1, . . . , n has action 

S ((/(■), 9(0) = / T £(g(t),g(t))dt (8) 
Jo 

and that this action is invariant if q % changes by 

5q\t) = e a (t)B?Mt)) (9) 

where for a — 1, . . . , m each e a is an arbitrary and independent function of 
time t and the summation convention has been used, so that 

= 8S 

e dtf{t) 
+ (e a WK(9(t)) + e a (t)^K(9(t))) ?^^ dt. (10) 

Then, since each e a is an arbitrary function of t, 

d£(q{t),q(t)) n 



K(q(t)> 



dq^t) 



and RMt)) wit) + dt RMt)) dm ~ 



(ii) 



(12) 



From (TlT!) we deduce that on taking the Legendre transformation there will 
be constraints 

(j) a (p, q) = R l a (q)Pi = a = l,...,m (13) 
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while f)12p shows that, when the position variables q l (t) satisfy the equations 
of motion 

d d£(q(t),q(t)) _ d£(q(t),q(t)) (u) 
dt 8qi(t) dq^t) ' 1 ' 

the Noether current J a = Rl(q(t)) ^—^ rrr- — is conserved, so that, if 

aWV dq l (t) 

(Pi(t) , q % (t)) is a trajectory in phase space which satisfies the equations of 
motion, then 

^foa(Pi(*),ff*(t)))=0- (15) 

Standard arguments show that the equations of motion ffT4l) for the La- 
grangian system may be regained from the equations of motion of the system 
with extended Hamiltonian 

H e (p, q, A) = H(p, q) + \y a { Pl , q { ) (16) 

where the variables X a ,a = 1, . . . , m are Lagrange multipliers for the con- 
straints and H is the canonical Hamiltonian 

H(p,q)=q i p i -£(q,q) . (17) 



In the Hamiltonian setting ( fl5l) corresponds to {H, <f) a } = . If additionally 
the constraints satisfy a relation of the form 

{^ a ,<p b } = C c ab <p c (18) 

then they correspond to the momentum map of an infinitesimal Hamiltonian 
action on the phase space R 2n = T*(IR n ) of a Lie group G whose Lie algebra 
referred to a particular basis takes the form ffl8|) . The constraints <p a generate 
infinitesimal transformations on phase space which map a solution of the 
equations of motion of the system to another solution. Such transformations 
are known as gauge transformations and make explicit the redundancy of the 
system which follows from the symmetries in the Lagrangian. Assuming that 
there is actually a Hamiltonian action of the full group G on the phase space, 
then the reduced phase space, which carries the true degrees of freedom of 
the system, is C/G where C is the submanifold of M. 2n where <j) a {p, <l) — 
0, a — 1, ... ,rn. A standard result of Marsden and Weinstein [7] shows 
that C/G has a natural symplectic structure. (The corresponding Poisson 
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brackets correspond to Dirac brackets if gauge- fixing functions exist.) In this 
ideal situation the BRST procedure together with appropriate gauge-fixing 
leads to a path integral quantization of the system in a manner which is well 
understood [U [H]. The BRST charge is 



where the ghosts if and ghost momenta 7r a are canonically conjugate anti- 
commuting coordinate functions on a super extension of the original phase 
space with super symplectic form typically 



(This construction can of course be expressed in intrinsic language.) The 
space of observables on the true phase space C/G can then be identified 
with the cohomology of Q at ghost number zero (where each ghost rf has 
ghost number 1 and each ghost momentum 7r a has ghost number —1). The 
quantization of the theory can be carried out using path integrals in the 
super phase space using the extended Hamiltonian H + [Q, x] where x is a 
gauge-fixing fermion for the theory. Path integration leads to a supertrace 
which can be shown to give a trace over BRST cohomology provided that x 
is chosen in such a way that the cohomologically trivial term [Q, x] regulates 
the infinite series in these traces [5]. 

While the procedure described above gives the basic principle of the BRST 
construction and its use to generate path integrals, many systems require an 
extension of these ideas. In the following sections the construction is extended 
to systems with secondary constraints, where it is necessary to consider the 
space of trajectories in phase space. 

3 Gauge generators and secondary constraints 

In this section systems whose symmetry leads to secondary constraints are 
considered. Starting from the Lagrangian picture, where the nature of sym- 
metries which lead to secondary constraints is identified, the Noether currents 
and primary constraints (which in this case are not the same) are identified, 
and then, largely following Pons [TU] , the canonical analysis of such systems is 
described, and a key condition f H6|) determining symmetries of the equations 



(19) 



dpi A dq 1 + dir a A dr] a . 



(20) 
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of motion is derived. Finally it is shown how solutions of this condition may 
be derived. This section draws on work by Castellani [5] and by Henneaux, 
Teitelboim and Zanelli [9j, as well as the work of Pons already cited. 

Suppose that as before we have a classical system with n degrees of freedom 
q\ i = 1, . . . , n and Lagrangian £(g, q), but that the action 

£(?(■).*(■))= f £(q(t),q(t))dt (21) 
Jo 

is invariant when the infinitesimal change in the variables q % takes the form 

Sq\t) = {e a (t)Ri(q(t)) + e a (t)^(g(t))) , (22) 

so that both arbitrary infinitesimal parameters e a and their time derivatives 
are involved. In this situation the equivalent of fflUj) is 

= 5S 

+ ^ a mi(m + e a (t)^ t K(q(t)) 

+ msMt)) + ?(t)±s<Mt)j) d£{ lflf )] • (23) 

Taking coefficients of e, e and e in turn we deduce that 

^))«f^ = ( 24) 

(*m) + ±S> Mt)) ) - 0(25) 

flnH Ri(nM) d£(q(t),q(t)) | d , » ^(g(t),g(t)) _ 

and + dt RMt)) ~ °" 



(26) 



From (J2"6]l we deduce that when the equations of motion (TTJJ hold 
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so that the Noether current i?^(g(t)) ^^M^ is conserved on shell as be- 
fore, while (1211) shows that on passing to phase space using the Legendre 
transformation there are primary constraints 

0« ee Si( q ) Pi = (28) 

but from ( |25l) we see that these are not preserved in time by the equations 
of motion, instead we have 

^(0i 1) (^),?W))=0i 2) (pW,g(t)) (29) 

(2) 

where (pa {p^q) = R l a (q)pi. The system thus has secondary constraints 
(j>a (p, <?) which, by ( l27|) . are preserved in time by the equations of motion, 
so that there are no tertiary or higher constraints. This analysis makes it 
clear that the signature of a system which leads to secondary constraints in 
the Hamiltonian setting is that the infinitesimal variations in the paths q l {t) 
which are symmetries of the action involve time derivatives of the infinites- 
imal parameter. (Extending the arguments above shows that if the action 
is invariant under changes in the fields which involve second derivatives of 
the parameters with respect to time, then there will be tertiary constraints, 
and so on. For simplicity in the remainder of this paper we will concentrate 
on systems with primary and secondary constraints but no higher order con- 
straints.) 

As before the equations of motion for the Lagrangian system (TH|) may be 
regained from the equations of motion of the system with Hamiltonian 

H e (p, q, A) = Vl q l ~ q 1 ) + \ a (j)£\p, q) • (30) 

(It is important to note that this extended Hamiltonian does not include the 
secondary constraints of the system, a point made very clearly by Pons [TO].) 
The equation of motion for X a (t) is of course 

<f>£\p(t),q(t)) = (31) 

so that the condition that the primary constraints be zero is included in the 
equations of motion. From f )28|) and f )29|) we deduce that 

{H e ,^}=^ and {tf e ,0i 2 )}=O. (32) 
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In the case of the toy model (pQ) the extended Hamiltonian takes the form 

H e (pi,P2,q 1 ,q 2 ,X) = \{P2) 2 -P2Q 1 + Xpi (33) 
leading to equations of motion 



P2 = 0, 
Pi = P2, 



e 


= P2 


e 


= A, 


pi 


= 



(34) 



which are equivalent to 

In many cases a covariant derivative is more appropriate, with the infinites- 
imal symmetry transformation of the path q % (t) expressed as 

5q\t) = e a (t)Rl(q(t)) + ^S'Mt)) (35) 

where 

Vt dt b\H\JJ 

for some connection M. In this situation analogous arguments to those 
leading to ([28]), (|27J) and (J29]) give 

£ ( Km fJMMM) = o (36) 



Vt V dtf(t) 

-(g(t),g 

&f(t) 



so that the Noether current R^(q(t))^S%$^- is covariantly conserved on 



shell, while as before we have the primary constraints </>i = S l a (q(t))pi = 
which satisfy 

^(^(PW, ?(*))) = tf } (p(*), <?(*)) (37) 
where (jy a = R\{q(t))pi. Additionally equations ( 152]) now have the form 

{^e,0i 1) }=0i 2) and {^,0i 2) } = -M a Vl 2) . (38) 

We will now investigate the gauge transformations which this system admits. 
Gauge transformations are mappings of the space of trajectories in phase 
space into itself which map each solution of the equations of motion with 
given initial conditions p l (0) = Pio,q , '(0) = q l for the variables Pi{t),q l (t) to 
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another satisfying the same initial conditions. Such solutions are not unique 
because of the presence of the Lagrange multipliers in the equations of mo- 
tion. (In the next section it is shown that such transformations are canonical 
in the symplectic manifold of paths in phase space, with infinitesimal gener- 
ators acting by Poisson bracket.) 

Now, largely following Pons [10] , we identify the condition which an infinites- 
imal generator G(t,pi(t), q l {t)) must satisfy if it is to generate a gauge trans- 
formation. Suppose that g(p, q) is an arbitrary function on phase space so 
that its evaluation on a phase space trajectory (p(t),q(t)) is g(p(t),q(t)). 
For any infinitesimal change of trajectory Pi{t) !->■ p%(t) + 5pi(t), q l (t) i-)- 
q l (t) + 5q l (t) we have 



Suppose now that (p;(t), X a (t)) and (pi(t) + 8pi(t), q\t) + Sq'ty), X a (t) + 
5\ a (t)) are both solutions of the equations of motion, and that 



so that G(t,p(t),q(t)) is the generator of the infinitesimal transformations. 



5^-g(p(t),q(t)) = ^-5g(p(t),q(t)) 



(39) 



Sg(p(t),q{t)) = {G,g}(t,p(t),q(t)) 



(40) 



Then 



%-g(p(t), q(t)) = {g, H e } (p(t), q(t), X(t)) 



(41) 



giving 



S ^9(p(t), q(t)) = {G, {g, H e }} (t,p(t), q(t)) + {g, 5H e } (p(t), q(t)) , (42) 
where 8H e {p{t),q{t)) = 5 \ a (t) ^ \p{t) , q(t)) . Also 
^-5g(t,P(t), q(t)) = {{G, g} , H e } (t,p(t), q(t)) + {^,g) (t,p(t), q(t)) (43) 




so that (using (1391 ) 




(44) 



Since g is an arbitrary function on phase space we must have 



— + {G,H e } + 5\ a <j ) P = f(t). 



(45) 
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Following Pons, we may absorb f(t) into G(t) so that finally we deduce that 
the requirement we can impose to determine a gauge generator is that 



dG 
dt 



+ {G 1 H e } + 5\yM = 0. 



(46) 



Also, so that the initial conditions of the trajectory are unchanged we im- 
pose G(0) = 0. The key condition (|46p is unchanged when the primary and 
secondary constraints are defined using the covariant derivatives of the in- 
finitesimal parameters. We now show how generators may be constructed 
which satisfy ( 146]) with this initial condition. (This condition will also deter- 
mine the change in the Lagrangian multipliers.) For simplicity (and because 
these seem to be the examples which appear in practice) we largely consider 
systems which require no more than secondary constraints, but we allow 
the more general situation of equation fl38|) which arises when the infinitesi- 
mal symmetry transformations of the Lagrangian involve the covariant time 
derivative of the infinitesimal parameters. The extension to tertiary and 
higher constraints does not present new difficulties. (The BRST symmetries 
will in fact apply to a yet more general situation.) 

For a system with primary first class constraints </>i^ (pi, q % ), a — 1, . . . , m and 
no secondary constraints, the result is well known: we set 



where each e a , a = 1, . . . , m is an arbitrary function of time with e a (0) = . 
In this case, 



— + {G,H e } + SX a (t)^ip i (t),q i (t)) 

= k\t)^\ Pl {t), q \t)) + e a {^\H e }+5\\t)^\ Pl {t), q \t)) (48) 



so that (H| is satisfied if 5\ a = -e a . 

When there are both primary constraints (p^\p,q),a = l,...,m and sec- 
ondary constraints (fr a (p, q), a = 1, . . . , m satisfying ( |38l) then we set 



G(t, Pl (t),q\t)) = e a (t)^\p(t),q(t)) - e a (t)^\p(t),q(t)) (49) 




(47) 
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with e a arbitrary function of time with e a (0) = 0. In this case 
BC 

— + {G, H e } (t,p(t),q(t)) + 5\ a (t)^\ Pl (t),q%t)) 

= ^{t)^\p(t),q{t))-e\t)^\p(t),t)) 

+e a (t) {^\H e } (p(t),q(t))-e(t) H e } + 5X a (t)^\ Pi (t), q\t)) 

(50) 

so that ( H6|) is satisfied if 5X a = — e. (A similar construction is possible in 
the covariant case.) 

For the toy model with Lagrangian (OQ), this prescription gives 

G = epi - ep 2 (51) 

so that, if 5 A = — e, 

BC 

— + {G,H e } + 5X Pl 
= epi - ep 2 + ep 2 + 8\pi 

= 0. (52) 

In the following section we show that by considering the natural symplectic 
structure on the space of paths in phase space we may recognise the gauge 
generators G(t) which are built from the primary and secondary constraints 
as moment maps. 



4 The symplectic geometry of paths in phase 
space 

In this section we briefly describe the symplectic structure which naturally 
exists on the space of paths on a symplectic manifold (following Wurzbacher 
[12]) and define a Hamiltonian group action on such a manifold together 
with the associated moment map. This allows us to recognise the gauge 
generators constructed in the previous section in terms of the moment map 
of an appropriate group action. 

Let M. be a 2n- dimensional symplectic manifold with symplectic form u, and 
local canonical coordinates pi, x l , i = 1, . . . , n so that 

u = dpi A dx { . (53) 
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Also let I = [0, T] be a non-empty interval in R and let At denote the set of 
smooth paths z : I — > At. This space can be given the structure of a Frechet 
manifold in a standard way, so that the tangent space T z At at a point z in 
Ai can be identified with the space C°°(I, z*TM). Thus a tangent vector X 
at z is identified with a vector field along z in Ai, so that in local coordinates 
X may be expressed as 



x = x'Mt)) A 



z(t) °Pi 



(54) 

*(t) 



It is then clear that there is a symplectic form uj on AA with 

Q(X, Y)(z(.)) = Xi(z(-))YP(z(.)) - Yt(z(.))X?(z(.)) , (55) 
so that the Poisson bracket of two functions F and G on At takes the form 

1 ' 1 5x*(z(-))5 Pi (z(-)) 6pi(z(-)) 8&{z{-)) • { ) 

A detailed account of symplectic structures on loop spaces may be found in 



It is now possible to define a symplectic action of a Lie group on At in 
the usual way, and also to define such an action to be Hamiltonian if the 
infinitesimal action is obtained by Poisson bracket with the transpose of the 
moment map and there is a well-defined symplectic quotient. Two simple 
examples will now be given, the second of which will show that the symplectic 
quotient of a Hamiltonian group action does not itself necessarily have the 
form of the space of paths in some symplectic manifold. It is precisely this 
more general situation which is required to handle second class constraints. 

The first example we give of a Hamiltonian group action on At arises when 
there is a Hamiltonian action of a Lie group G on the symplectic manifold 
At with momentum map J. If G denotes the group of paths g : I — >■ G, with 
group composition defined pointwise, then the action of G on At defined by 

g : z H- zg with zg(t) = z{t)g{t) (57) 

is readily seen to be Hamiltonian with moment map 

J:M^g* with J(z(-))(t) = J(z(t)) . (58) 
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Here the dual g* of the space g of paths in the Lie algebra g of G is defined 
to be the set of smooth paths a : I — > g*, the action of such a path on a 
smooth path rj : / — > g being defined by 

(a,ri)= [ dt(a(t), V (t)) . (59) 
Jo 

In this case if B is the symplectic quotient of Ai under the Hamiltonian 
group action of G then the symplectic quotient of Ai is simply the space B 
of paths in B. 

A second example of a Hamiltonian group action emerges from the toy exam- 
ple fll]). Here Ai = M 4 with canonical coordinates (pi,p 2 , x 1 , x 2 ) and G — R, 
the group of paths in the additive group R. Denoting a typical element 
e : / — > M. of 1R by e(-), the action of R on M. is defined by 

e(-) : ( Pl (-),p 2 (-),x\-),x 2 (-)) 

^(p 1 (.),p 2 (-)^ 1 (-) + e(-)^ 2 (-)-e(-))- (60) 

The Lie algebra g* may also be identified with the space of paths in R, and 
the transpose T of the corresponding moment map is defined by 

T(e)(t)=e(t) Pl (t)-e(t)p 2 (t). (61) 

In this case the symplectic quotient does not take the form of the space of 
paths in a finite-dimensional symplectic manifold. 



5 BRST cohomology and secondary constraints 

In order to develop a BRST formalism for the more general symmetries 
which lead to secondary and higher constraints an alternative, equivalent 
approach to phase space and the corresponding BRST operator is required; 
this approach will first be described for a system with no symmetry and then 
extended to the general case. This approach automatically generates, and 
justifies, the ghost terms which must be added to the Hamiltonian when a 
system has secondary constraints. 

Suppose that H is the Hamiltonian of a dynamical system on the symplectic 
manifold Ai which has no symmetries (and hence no constraints). Then 
consider the symplectic manifold hA = Ai x M? with symplectic form 

u = u + dE A dt (62) 
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where E and t are global coordinates on M 2 Woodhouse Now if we 

impose on M. the constraint <$e = H — E = Owe are in fact using the 
constraint which corresponds to the moment map of the following action of 
the additive group K. on M.: 

RxM^M, (t,(s,z))i-> (s + t,p?(z)) (63) 

where pf{z) is the flow of the Hamiltonian vector field Xh = fp-^r — 
starting from z, so that pf ' (z) is the solution at time t of the equation of 
motion with initial conditions z. 

Now, following the standard Marsden Weinstein reduction process, we see 
that if C is the submanifold of A4 where H — E = 0, then C/R is the space 
of trajectories z : R — > Ai which satisfy the equations of motion for the 
Hamiltonian H. This can of course be identified with A4, given that there is 
a unique solution of the equations of motion starting at each point of A4. 

We will now construct the BRST cohomology for this constrained system 
and show that it does lead to the quantum dynamics of the system. The 
standard construction for this single constraint gives as BRST operator 

= G (H - E) (64) 

where a is the unique ghost corresponding to the basis { ^ } of the Lie algebra 
of the additive group K. If we now choose as gauge-fixing fermion \ = rt 
where r is the canonical dual to o (so that the symplectic form on the super 
phase space is uj + dr A da) the gauge-fixing term is 

P,X) = (H-E)t + T<T. (65) 

On quantization the term — Et transforms from the Heisenberg picture (in 
which observables carry the dependence on time ) to the Schrodinger picture 
with states tp(q,t,r), = —i-^j^a = Path integration in the standard 
way will evaluate the supertrace of the operator exp (— Ht + rcr); direct cal- 
culation shows that the supertrace of exp (rcr) is simply 1 (as is also the case 
for exp (rcr + Aa) provided that A is independent of r). As a result carrying 
out cr and r integrals in the path integral will simply give a factor 1, and 
the remaining path integral will simply be the standard path integral for the 
trace of the imaginary time evolution operator exp (—Ht), as of course would 
be expected. 
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At this stage we have regained a standard result by an apparently pointlessly 
roundabout route; it will now be shown that this approach can be extended 
to handle systems whose symmetries lead to secondary and higher order 
constraints, deriving the extended Hamiltonian including ghost terms. 

Suppose that we consider the system whose classical symmetries were anal- 
ysed in Section [3j Recall that this system has Hamiltonian H(p,q), primary 
constraints <pa and secondary constraints <p^a \ with infinitesimal generators 
of classical symmetries taking the form G(t) = e a (t)<j)^ + e(t) a </>i 2 ' ) , corre- 
sponding to a group action of G on Ai. This group acts on the trajectories 
of the Hamiltonian vector field Xjj, and so we can consider the correspond- 
ing BRST procedure. Our parametrisation of the space of these trajectories 
as the quotient space of all paths in M. under the action of R generated by 
H — E means that we identify trajectories by their position at some arbitrary 
but fixed time t . Because at the particular time t the parameters e a (t ) and 
e a (t ) are all independent we effectively add the constraints </4^ and <j)^ each 
with its own ghosts and can construct the BRST operator in the standard 
way. At this stage it seems simpler to abandon a notation which distinguishes 
between primary and secondary constraints and simply consider constraints 
ip a and relations 

{V^, M = C c ab ^ c , {H, ^} = Vfy b (66) 

where each C c ah is a constant but V£ will in general depend on p, q. The 
corresponding quantum BRST operator is then 

ft = n H + ft G (67) 

with ft^ = a{H{p, q)-E) + ar] a 7c b V^ and ft 9 = r/>a + |Q^VVv (When V* 
is not constant a slightly more complicated BRST operator may be required.) 

The BRST quantization is completed by choosing gauge fixing fermion \ = 
Xh + Xg where xh = i~t and xg is chosen to have the necessary features of a 
gauge-fixing fermion for ft^ [5]. Using the Et term as before to move to the 
Schrodinger picture, we obtain the quantum Hamiltonian 

p,x\=Ht + V^ a n b t + aT + Acr+ pa, xg] (68) 

where A is independent of r. Once again we can integrate out a and r, so 
that the remaining path integral will evaluate the supertrace of the evolution 
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operator exp [H t + V^rf-K h t + [Q G , Xg})- This is a standard result, but the 
arguments we have given fully justify the addition of terms V^rj a 7ib to the 
classical Hamiltonian, rather than simply including them in order to obtain 
a Hamiltonian which commutes with Qq- 

To end this section this procedure will be carried out for the toy model with 
Lagrangian ([T|). Using (I5TI) and (|67|) we see that in terms of two ghosts rj 1 , if 
the BRST charge is 

Vt = n G + tt H (69) 

where 

tt H = a(H -E) + ari\ 2 (70) 
= rfpi + r] 2 p 2 . (71) 

In order to construct a gauge-fixing term for this theory we bring in non- 
minimal gauge fields and symmetries as first introduced by Fradkin and 
Vilkovisky [IB]. These fields can be interpreted as giving dynamics to the 
Lagrange multiplier, or alternatively simply introducing extra degrees of free- 
dom which are then eliminated by a further symmetry. When no secondary 
constraints are involved a new field l a with conjugate momentum k a is in- 
troduced for each constraint <f) a , while the Lagrangian is unchanged, so that 
each k a is constrained to be zero. The corresponding BRST operator then 
has additional terms n a k a where K a is the ghost corresponding to the con- 
straint k a = 0. In the case of secondary constraints extra terms must be 
added to the Hamiltonian. For the toy model there are two constraints and 
so two new fields I 1 , 1 2 with conjugate momenta k ly k 2 are introduced and 
an extra term -l l k 2 is added to the Hamiltonian, so that k\ is a primary 
constraint and k 2 a secondary constraint. (Here e is a constant factor.) The 
corresponding BRST charge is then 

Ojvm = Qh + Qg + ^k where Qk = +-<yl 1 k 2 + -an 1 X 2 + K 1 ki + K 2 k 2 (72) 

e e 

and Ai and A2 are the ghost momenta conjugate to k 2 and k 2 resepectively. 
It will now be shown that choosing the gauge- fixing term 

Xnm = rt + nil 2 + -q Y \ 2 (73) 
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leading to the gauge-fixed Hamiltonian 



H NM = H c + -l l k 2 + 77V2 + 



^2 

e e 



+ttiK 2 + -A^ 1 + -q x k 2 + / 2 pi (74) 

after the Et term has been used as before and integration over a and r 
has been carried out. Inspecting this Hamiltonian we see that it has the 
standard form to ensure that the supertrace of exp (iH NM t) will project to 
the trace over physical states [5]. This may be seen from the fact that 
H NM = HD + lml 2 + ig 1 A 2 ,n G + ^] where H D = H c + \l 1 k 2 + r} 1 ir 2 + \n l \ 2 
and by construction [Hp, + ^k] = together with the fact that {k\1 2 + 
^q 1 X 2 , &>g + &>k] — + ^-A^ 1 + \q x k 2 + l 2 pi which is the standard form 
for the quartet mechanism [3J. 

For completeness, it will now be shown that the corresponding path integral 
reduces to the form given in [13] after integration over some of the variables. 
The path integral takes the form 



Vp 1 Vq 1 Vp 2 Vq 2 Vk l Vl 1 Vk 2 Vl 2 ViT 1 Vr] 1 V7r 2 Vr] 2 VK 1 V\ l VK 2 VX 2 
exp^y i{piq l +p 2 q 2 + kj 1 + k 2 l 2 + iiif] 1 + n 2 rj 2 + « 1 Ai + k 2 \ 2 ) 

+H NM )dt] . (75) 



To evaluate this path integral we use a technique first introduced by Fradkin 
and Vilkovisky [16] and scale k 2 — > ek 2 and A 2 — > eA 2 (which leaves the 
measure unchanged) and let e — > 0. The integrand then becomes 

exp^y i{piq l +p 2 q 2 + kj 1 + nxrj 1 + 7r 2 ?) 2 + k^Ai) 

+~J?l - Q l P2 + 1% + V^2 + + vtxk 2 + A^ 1 + g 1 ^ + i 2 pi)dij (76) 

Integrating out Z 1 gives a factor 5(/c 2 — ki), while integrating out I 2 gives a 
factor <5(pi). On the ghost side, integrating out 7r 2 , n l and k 2 give factors 
5{jf + 77 1 ), 5(X 2 — Ai) and (^(tt 1 ) respectively. Integrating out pi,p 2 , k 2 , tti, 77 1 
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and A2 is then straightforward, reducing the path integral to the form 
Z = J Vq x Vq*Vk{Di?V\ x 

expQTi Qte 1 - ?? + hq 1 + dtj . (77) 

in agreement with (4.16) of [13] if k\ is replaced by b, Ai by C and rj 2 by C. 

6 Examples 

In this section we present some examples to illustrate how primary and sec- 
ondary constraints are combined in a single generators of gauge transforma- 
tions that satisfies Pons' condition (jl6"j) . These examples are the relativistic 
free particle, pure Yang-Mills theory and the bosonic string, all of which have 
secondary constraints. The first of these is an abelian theory, whereas the 
second and third are non-abelian field theories. For each example, the action 
is presented. Then, following Dirac's algorithm, the constraints are deduced 
and a generator of the gauge transformations is constructed which is shown 
to satisfy the required condition. 



6.1 Relativistic free particle 

This system corresponds to a (0+l)-dimensional gravity theory with scalar 
fields. The degrees of freedom are particle coordinates and an einbein e, 
each of which is a function of a single proper time variable r. The action for 
this system is given by 

S = |dr£=i/drp4l^-m> e (T)) , (78) 

where a dot over a variable indicates a derivative with respect to proper time 
t. From the action (1T8]) . the conjugate momenta to e and x M are respectively 

P = — = and P = — = ^ (79) 

so that we have the primary constraint 

0(1) = p e = . (80) 
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The extended hamiltonian is then 

H e = H + v 1 cj )l , (81) 

where v 1 is an arbitrary function of r, and H is the canonical Hamiltonian 
which takes the form 

H = l -e{P^ + m 2 ) . (82) 
Since = § (P M P" + m 2 ) we have the secondary constraint 

0( 2 ) = ^(P^ + m 2 )=O . (83) 

This new constraint satisfies {if, f^ 2 **} = and so there are no further con- 
straints. 

Finally we construct the gauge generator G(t) = ecf)^ + ecj)^ which satisfies 
()46l) provided that 5X = e. 



6.2 Yang-Mills 

This is a very well known gauge theory where the fundamental fields are 
gauge potentials A^,a = 1, . . . , dim G where G is the gauge group of the 
theory. The Lagrangian density for this system is 

C = \f* v F^ (84) 

where F* = d^A® — d u A^ — gf% c A b A c v with g the coupling constant and 
f a bc the structure constants of the Lie algebra of the gauge group G. The 
conjugate momenta are 

<-§r° (85) 

< = H = -n ( 86 , 

so that there are primary constraints 

= < (87) 
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while the canonical Hamiltonian is 
H = |d 3 x Qtt'Cx.Ot&x.O-^ 

(88) 

where the covariant derivative is defined by T> il r K v a = d^n^ + gA b f^tr" ■ and 
spacetime has been taken to have four dimensions. 

Since {H,<j>^'\ = T>iir l a we have secondary constraints 

€ ] = -A< = (89) 

These new constraints satisfy 

{H,^} = -gf abc A b )( j ) 2 c ^0 (90) 

and so there are no tertiary or higher constraints. The Poisson brackets of 
the constraints with each other are 

{$\x,t),$\y,t)} = 0, {0«(x,t),0[ 2) (y,t)} = O 

and {0i 2) (x,t),0i 2) (y,t)} = / a c b 0i 2) (y,^(x-y). (91) 

The form of the constraints suggests the generator 

G(e) = (V) <Pl + e> 2 (92) 

with each e a an arbitrary function of time. Now 

d G = (<9 V - gf abc A c d e b ) <\>\ + (d e a ) <f> 2 a (93) 
and {G,U C } = -d ey 2 a (94) 
so that 8 G + {G, H c ] = (V d e a ) <\>\ (95) 

and thus (j46H is satisfied if 5X a = Po<%e a . Direct calculation shows that 

{G{e 1 ),G{e 2 )} = G{e 3 ) (96) 

where = d\t\f bc so that the symmetry group is the space of maps of 
spacetime into G, but the group action is not pointwise. 
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6.3 Bosonic String 

Relativistic free strings propagating in an arbitrary D- dimensional fixed back- 
ground space-time can be described with the Polyakov action [T7| 

5[ 7 ab ,X^] = a J d^^r'daX^X'g.uiX) (97) 

where X M are back-ground spacetime coordinates, g^ v is the background 
spacetime metric, 7 ab is the world sheet metric and £ stands for the two 
worldsheet coordinates r and a. When r y ah is taken in the ADM form [TB] 



i ab = ( r,_^i (98) 

V N 2 x N 2 J 

the Lagrangian density in the Polyakov action takes the form [19] 

C = f-X^g^ + 2\X»X"'g^ + ( Y " A ') X '" X ' V 9^ (") 

where 

dX^ dX^ 
The conjugate momenta to X^ are: 

p^-^j^+^jpv (ioi) 

while those of X,N,x are 

P x = P N = P x = (102) 
giving the primary constraints 

<p [ x ] = Px, <P ( n ] = Pn and ^ = P X . (103) 
The canonical Hamiltonian is 

H c (t) = J da H(x(a,T),N{a,T),X(a,r),X^a,T),P^a,T)) (104) 
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where 



4 2) 



direct calculation shows that 1 19 



where k = e'a — ole. 
Since 



(105) 



Taking Poisson brackets with the primary constraints gives 

{H C ,^} = -X^P, (106) 
{H c , <^ } } = -^/= {P,Pu9^ + We 2 X'»X'»g, v ) (107) 

{He = ~^j- 2 {P^r + 4a 2 e 2 X'-X'^) (108) 
and hence we have secondary constraints 



X'*Pp (109) 
<f>P = P„P v gT + laWx'^g^ . (110) 

If we define 



jdae(a)^\a) and $\e) = J dae{a)<j>%\a) (111) 



{A 2) (4< 




:«)> = 


A 2) («) 


(112) 


{A 2) (4< 




:«)} = 




(113) 


{<t>%\e),< 




»} = 


0g } ( K ) 


(114) 



H c = -p^P + A0g } (115) 

(where p = 4a gV^ ) there are no tertiary or higher constraints, which suggests 
a generator of the form 



+e D (a,T)^\a,T)+e H ^\a,r)) (116) 
24 



with each e an arbitrary function of r and a. By substituting this generator 
in the equation (EES]) , we find 

dG 



dr 



{G,H} 



da 



'TV 



+ (e° - + (A'e° - A(c°)') - (p'e H - p(e")')) 0g } 



! ( e" + ^ - ^ - {p>e» - p(e»)') + (AV + A(e")' 



so that the parameters must obey the relationships 



t(2) 



117) 



+ (AV - A(^)0 + (p'e* - p(e 



(118) 

t N = f*-- (e H - (A'e* + X(e H )') + (p'e D - p{e D )')) (119) 



and the generator is 



G 



da 



(e D + (A'c - A(c D ) / ) + (p'e H - ^ 
+ {^f ~^ H - (A'e* + A(e*)') + (p'e^ - p(e°)'))) 
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